■’^’Calhoun 

IflitiKuliarttl  Arehivtaf  (he  Naval  ftutgrdduH*  School 


Authors) 

Overbay,  WilliamA. 

Title 

Sturmfs  theorem  and  the  zeros  of  a  sol  uti  on  to  a  different]  al  equati  on. 

Publisher 

Monterey,  California:  U.S.  Naval  Postgraduate  School 

Issue  Date 

1963 

URL 

http://hdl .  handl  e  net/10945/13209 

This  document  was  downloaded  on  May  12,  2015  at  03:37:29 


DUDLEY 

KNOX 

LIBRARY 


htt  p:/ /www.  n  ps.  e  d  u/l  ib  ra  ry 


Calhoun  is  a  project  of  the  Dudley  Knox  library  at  NP5,  furthering  the  precepts  and 
goals  of  open  government:  and  government  transparency.  All  information  contained 
herein  has  been  approved  for  release  by  the  NPS  Public  Affairs  Officer. 

Dudley  Knox  Library  /  Naval  Postgraduate  School 
411  Dyer  Road  /  1  University  Circle 
Monterey,  California  USA  93943 


http://www.n  ps .  edu/ 


NPS  ARCHIVE 
1963 

OVERBAY,  W. 


STURM’S  THEOREM  AND  THE  ZEROS 
Of-  A  SOLUTION  TO  A  DIFFERENTIAL  EQUATION 

WILLIAM  A.  OVERBAY. 


STURMS  THEOREM  AND  THE  ZEROS  OF  A  SOLUTION 
TO  A  DIFFERENTIAL  EQUATION 


*  *  *  *  *  * 


William  A0  Overbay 


STURMS 


THEOREM  AND  THE  ZEROS  OF  A  SOLUTION 
TO  A  DIFFERENTIAL  EQUATION 


by 

William  A.  Overbay 

// 

Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 


with  major  in  mathematics 

United  States  Naval  Postgraduate  School 
Monterey,  California 


library 

u.s.  naval  postgraduate  school 


STURMS 

MQN  f  ;*pcy  CAI.JPORNIA 

THEOREM  AND  THE  ZEROS  OF  A  SOLUTION 

TO  A  DIFFERENTIAL  EQUATION 

by 

William  A.  Overbay 

This  work  is  accepted  as  fulfilling 

the  thesis  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

with 

major  in  mathematics 

from  the 

United  States  Naval  Postgraduate  School 


AB  ST  RAC  T 


The  number  of  zeros  contained  in  a  given  interval  (a,b)  for  a 
solution  to  a  Sturm- Liouvi lie  differential  equation  is  of  importance  in 
many  problems  of  mathematical  physics.  This  number  maybe  deter¬ 
mined  through  Sturm* s  Comparison  Theorem.  Given  one  zero  of  the 
solution  to  a  Sturm- Liouville  differential  equation,  a  technique,  based 
upon  Sturm*s  Theorem,  of  computing  the  next  consecutive  zero  of  the 
solution  is  proposed.  The  existence  of  a  function  which  satisfies  the 
desired  end  results  of  the  proposed  technique  is  shown.  The  technique 
is  then  applied  to  Bessel’s  differential  equation  and  the  results  tabula¬ 
ted  for  the  first  20  roots  of  JQ(x)  and  Jj(x).  Unfortunately  this 
technique  did  not  achieve  the  desired  result  of  convergence  to  succes¬ 
sive  zeros  of  the  given  Bessel  Function. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance 
and  encouragement  given  him  by  Professor  E.  J.  Stewart  of  the  U.  S. 
Naval  Postgraduate  School  in  this  investigation. 
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1.  Introduction. 


Consider  a  differential  equation  system. 

[  r(x)  y'  (x)]  +  P(x)  y(x)  =  0  (1.  1) 

y(a)  =  0  y(b)  =  0 

where  r(x)  and  p(x)  are  positive,  continuous,  real  valued  functions 
on  [a,b].  The  system  (1.  1)  is  known  as  a  Sturm- Liouvilfe System. 

A  frequent  problem  encountered  in  mathematical  physics  requires 
the  determination  of  consecutive  characteristic  roots  of  a  system  (1.  1)l 
The  solutions  to  (1.  1)  usually  do  not  possess  the  nice  regular  charac¬ 
teristics  of  the  trigonometric  functions;  consequently  the  former  are 
in  general  more  difficult  to  manipulate. 

When  the  value  of  one  characteristic  root  a  of  (1.  1)  is  known, 
Sturm's  Comparison  Theorem  suggests  a  possible  means  of  computing 
the  next  consecutive  root  b  in  the  following  sense: 

A  sequence  of  functions  is  generated  from  (1.  1).  Each  element  of 
the  sequence  is  obtained  by  selecting  a  particular  value  of  the  indepen¬ 
dent  variable  x  =  x^  and  substituting  this  value  into  the  coefficient  P(x) 
and  r(x)  of  (1.1)  giving 

[r(x.)y'(x)]'  +  P(x.)  y  (x)  =  0  (1.2). 

The  solution  of  (1.  2)  is  forced  to  vanish  at  y  =  a,  and  the  next  consec¬ 
utive  zero  of  this  solution  is  chosen  as  the  next  element  of  the 

sequence.  It  is  desired  that  the  sequence  so  generated  result  in  an 
equation  of  the  form  (1.  2)  such  that  its  solution  vanishes  at  x  =  a 
and  x  =  b,  the  zeros  of  the  solution  of  (1.  1). 

The  theory  is  developed  for  a  general  Sturm- Liouvilfe  System  but 
specific  applications  are  limited  to  Bessel's  Equation  of  order  n  . 

(n  a  positive  integer). 
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We  shall  define  an  indexing  of  the  roots  of  a  solution  to  a  differen¬ 


tial  equation  as  follows: 

Let  rj  be  the  first  root  >  0,  and  r^  the  next  root  greater  than 
r  y  Furthermore,  if  r^  is  the  i  th  root  y  0,  then  r^+^  is  the 
(i+l)st  root  >  0. 

If  r^  is  the  given  root  of  a  solution  to  a  differential  equation,  we 

shall  define  r,  ,  ,  as  the  next  consecutive  root  of  the  solution. 
k+1 
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2.  Mathematical  Background. 


This  section  contains  theorems  and  mathematical  developments 

which  will  be  called  upon  in  later  sections. 

2,  1.  Sturm's  Comparison  Theorem. 

Given  two  equations  of  a  Sturm- Lao uville  type: 

y"+Ay=0  (2.1.1) 

y(a)  =  0 

u"  +  Bu  =  0  (2.  1.  2) 

u(a)  =  0 

where  A  and  B  are  positive,  real  constants  and  B  is  greater  than  A. 

Let  us  determine  a  solution  y(x)  to  (2.  1.  1)  such  that  y(a)  =  0. 

The  elementary  theory  of  differential  equations  tells  us  that  the  general 
solution  is  of  the  form 

y(x)  =  sin  '/Ax  +  C-^  cos  A  x  . 

When  x  =  a 


y(a)  =  Cj  sin  V^a  +  C2  cos-y^a  =  0 
C  i  sin'll  a 


C2  = 


cos 


yr. 


Since  is  an  arbitrary  constant,  choose  =  kQ  cos  VTa. 

y(x)  =  kQ  [sin  yfK  x  cos  Via  -  c  os  VTX  sin  VTa] 
y(x)  =  kQ  sin  'sj~A  (x  -  a). 

Similarly,  (2.  1.  2)  will  have  a  solution  u(x)  =  k^  sin  Vb(x  -  a)  . 

The  assumption  A  <  B  tells  us  u(x)  will  "oscillate"  more  rapidly  than 
y(x)  or  that  u(x)  will  have  at  least  as  many  zeros  as  y(x)  on  a  given 
•  ..interval  (a,  b). 

When  the  coefficients  of  (2.  1.  1)  are  positive,  real  valued  functions 

of  the  independent  variable  x,  the  solutions  obtained  must  be  defined 
in  terms  of  so-called  "special"  functions  such  as  Bessel  functions. 
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Charles  Sturm  (1803-1855)  made  a  study  of  the  rates  of  "oscilla¬ 
tion"  of  solutions  to  equation# of  the  form 
[  r  (x)  y'(x)  ]'  +  P  (x)  y(x)  =  0 
[  r  (x)  w'(x)  ]'  +  Q(x)  w(x)  =  0 
y(a)  =  w(a)  =  y(b)  =  Q 

where  r(x)^>0;  P(x)  and  Q(x)  are  continuous  on  [a,b].  He  proved 
that,  in  general,  the  larger  P(x)  the  faster  y(x)  will  oscillate.  More 
concisely,  Sturms  Comparison  Theorem  is  stated  as  follows: 

Theorem  2.1  Sturm's  Comparison  Theorem. 

Given  two  equations; 

[  r(x)  y'  ]  '  +  p(x)  y  =  0  (2.  1.  3) 

and  [  r(x)  v'  ]'  +  q(x)  v  =  0  (2.  1. 4) 

where  r  (x)  is  positive  in  the  closed  interval  [a,b]; 
r(x),  p(x)  and  q(x)  are  continuous  in  [a,b]  and 

q(x)  >  p(x)  with  strict  inequality  for  at  least  one  point  in 

[a,b].  If  a  solution  y(x)  of  (2.  1.  3)  has  consecutive 

zeros  at  x  =  x  and  x  =  x,  ,  with  x,  >  x  ,  and  x,  ,  x 
o  1  l^o  1  o 

are  in  the  interval  [a,b],  then  a  solution  v(x)  of  (2.  1.4), 
which  vanishes  at  x  =  xq,  will  vanish  again  in  the  open 
interval  (xQ,  x^)  . 

Proof: 

Suppose  y(x)  )  0  in  (xq,  Xj)  and  y'(xQ)>0,  y'(xj)  <  0  and 
v'(xo)>0.  Then  if  y(x)  and  v(x)  are  solutions  to  (2.  1.3)  and 
(2.  1.4)  respectively,  we  have  the  identities 

[  r(x)  y'(x)  ]'  +  p(x)  y(x)  =  0  (2.  1.  5) 

and  [  r(x)  v'(x)  ]'  +  q(x)  v(x)  =  0  (2.1.6) 
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If  we  multiply  equation  (2.  1.  5)  by  -v(x)  and  equation  (2.  1.6)  by 
y(x)  we  have 

■»  [r(x)  y'(x)  ]'  v  (x)  —  p(x)  y(x)  v(x)  =  0  (2.  1.  7) 

[ r(x)  v'(x)]'  y(x)  +  q(x)  v(x)  y(x)  =  0  (2.  1.8) 

Adding  (2.  1.  7)  and  (2.  1.  8)  then  integrating  from  xq  to  x^  we 
have 


,X1  ,-xl 


’(x)[ y(x)v'(x)  -  y'(x)v(x)]  I  +T  [q(x)-p(x)]y(x)v(x)dx  =  0 

o  J-x. 


or 


TX1 

rfx^y'UjMxj)  =  /  [q(x)-p(x)]  y(x)  v(x)  dx  (2.  1. 

/x 


9) 


o 

since  y(xQ)  =  0  if  xq  is  a  solution  to  (2.  1.3). 

Suppose  now  that  v(  x  )  >  0  everywhere  in  (x^,  x^),  (i.  e.  if 
v(xQ)  =  0  and  v'  (xQ)  >  0  then  v(x)  has  not  vanished  anywhere  in 
(xQ,  Xj)  ),  then  this  implies  that  r(x^)  v(x^)  y'(xj)  is  negative  while 


x 

(q(x)  -  p(x)  )  y(x)  v(x)  dx  is  positive.  This  contradiction  implies 
o 

that  v(x)  cannot  be  positive  everywhere  in  (xQ,  x^).  Similarly*  the 
assumption  v'(xQ)  <  0  and  v(x)  negative  everywhere  on  (a,b)  again 
leads  to  a  contradiction  and  the  theorem  is  established. 

2.2  The  Normal  Form. 

Consider  a  differential  equation 

y"  +  P(x)  y'  +  Q(x)  y=  0  (2.  2.  1). 

Let  y=  u  v,  then  y'  s  uv'  +u'v 


and 


y’fcuv"  +  2u'v'  +  u"v. 

Substituting  for  y  in  (2.2.1)  we  have 

uv"  +  [  2u'  +  P(x)u]  v'  +  [u»  +  P(x)  u'  +  Q(x)u]  v  =  0  (2.  2.  2). 
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Suppose  we  choose  u  such  that  the  coefficient  of  v1  vanishes.  Set 


2u*  +  P(x)  u  =  0 

du/u  +  ^2/^  dx  =  0 

u  =  C  exp  [-  dx] 

Let  C  =  1  and  we  have 
u*  =  j  P(x)  u 

u11  =  2’P(x)  u1  -  ^  u  ^^dx 

L  (P(x))2  1  dP/j 

=  +  — ~4~~  u  -  2  '  dx  u  . 

Substituting  u*  and  u"  in  (2.2.2)  we  get 


uv 


"  +  t-P(x)  u  +  P(x)  u]  v'  +  [Q(x)u  +  u  -  i  ^^*/dx  u 


P(x)‘ 


u  ]  v  =  0  . 


We  see  that  we  have  a  factor  of  u 
divided  out  leaving 


in  every  term, therefore  it  may  be 


v"  +  [ Q(x)  -  ^ - £  ]  v  =  0  (2.  2.  3) 

We  say  that  (2.  2.  3)  is  equation  (2.  2.  1)  transformed  into  normal 
form.  We  will  see  that  the  normal  form  plays  an  important  role  in  our 
later  development. 

2.  3.  Existence  Theorem. 

Consider  a  differential  equation  reduced  to  normal  form 


y"  +  R(t)  y  =  0.  (2.3.1) 

y(a)  =  y(b)  =  0. 

We  desire  to  know  if  there  exists  a  t  in  the  interval  (a,b)  such 
that  the  equation 

u"  +  R(F)  y  =  0  (2.  3.2) 

R(f)  is  R(t)  evaluated  at  t  =  F  such  that  a  solution  u(t)  to  (2.  3.  2) 
can  be  forced  to  have  zeros  at  both  t  =  a  and  t  =  b.  The  answer  to 
this  question  may  be  stated  in  the  form  of  a  theorem  as  follows: 
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Theorem  2.  2. 


Given  a  differential  equation  of  the  form: 

Y"  +  R(t)  y  =  0  (2.  3.3) 

where  R(t)  is  positive,  real  valued,  continuous  and  monotone  in  (a,c), 
where  b  <  c,  and  whose  solution  y(t)  has  consecutive  zeros  at  t  =  a 
and  t  =  b,  a  <  b  <  c.  Then  there  exists  a  t  in  the  interval  (a,b)  such 
that  the  equation 

u  *  *  +  R(F)  u  =  0  (2.  3.4) 

has  a  solution  u(t)  with  a  zero  at  t  =  a  and  moreover,  u(t)  has  its 
next  consecutive  zero  at  t  =  b. 

Proof: 

First  assume  R(t)  is  monotone  decreasing,  then 

R(a)  £  R(t)  >  R(b)  for  all  t  in  (a, b).  Choose  that  solution  u^(t)  of 

uM  +  R(b)  u  =  0  such  that  u^(t)  vanishes  at  t  =  a.  A  solution  that 

satisfies  this  requirement  is  u^ft)  =  sin  R(b)  (t-a).  Clearly 

u^(a)  =  0  and  moreover,  u^(t)  will  have  its  next  consecutive  zero, 

for  t  greater  than  a,  when  the  argument  of  the  sine  function  is  equal  tt, 

that  is ,  when  VSb)  (t-a)  =  tt  .  Since  R  is  assumed  decreasing, 

R(b)  <  R(t)  for  all  t  in  (a, b[ therefore,  Sturm's  Comparison  Theorem 

tells  us  that  u^(x)  will  vanish  again  at  t  =  t^  where  b  <  ty  Similarly, 

let  uM  +  R(a)  u  =  0  have  a  solution  u  (x)  which  vanishes  for  t  =  a. 

a 

Then  with  R(t)  <  R(a)  for  all  t  in  the  interval  (a,b),  have  that 
ua(t)  will  vanish  again  at  some  point  t  =  t^  such  that  a  <  t^  <  b. 

Case  II. 

Now  suppose  R(t)  is  increasing.  Choose  u^t)  a  solution  of 
v"  +  R(b)  v  =  0 

such  that  a  root  at  t  =  a.  R(t)  <  R(b)  for  all  t  in  (a,b),  hence  vfc(t) 
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will  vanish  again  when  t  =  t^,  <  b.  Similarly,  v"  +  R(a)  v  =  0  has  a 

solution  v  (t)  which  vanishes  at  t  =  a  and  will  have  its  next  consecutive 
a 

zero  at  some  t  =  tA,  b  <  t,. 

4  4 

In  either  case,  R(t)  increasing  or  R(t)  decreasing,  the  root  b  of 
equation  (2.  3.  3)  is  between  two  roots  of  (2.  3.4)  determined  by  evalua¬ 
ting  R(t)  at  t  =  a  and  at  t  =  b.  This  and  the  continuity  of  R(t)  assures 
us  that  there  exists  a  t  in  (a,b)  such  that  u^  (t),  a  solution  to  (2.  3.4) 
which  vanishes  at  t  =  a,  will  vanish  again  when  t  =  b,  and  the  theorem 
is  established. 

An  interesting  conclusion  from  Theorem  2.  2  is  the  following: 
Theorem  2.2.  1. 

Given  a  system  of  the  type  (2.  3.  1)  where  R(t)  is  monotone,  then 
the  next  consecutive  root  of  a  solution  u(t)  to  (2.  3.4)  such  that 
u(a)  =  0  must  lie  in  the  interval  (a  +  — - -  ,  a  =  -------  )  . 

V  R(b)  YR(a) 

The  proof  follows  directly  from  the  proof  of  Theorem  2.  2.  The  end 
points  of  the  interval  (a  +  "/yiUb),  a  +  U/  TitoT)  )  are  precisely  the 
t.  found  in  Theorem  2.  2. 

l 

Note  that  the  arbitrary  constants  occuring  in  the  solutions  to  equations 
(2.  3.  1)  and  (2.  3.4)  have  no  effect  upon  the  zeros  and,  therefore,  will 
be  ignored  in  all  further  considerations. 

We  should  also  note  that  in  a  general  Sturm- Li ouvilfe  System  (2.  3.  1), 
the  function  R(t)  is  not  always  monotone. 


y^+xy1  +  (x^-  n^)  y  =  0 


3.  Bessel's  Equation. 

Let  us  now  consider  a  typical  Sturm  -Liouville  System: 

(3.  1) 

y(a)  -  y(b)  =  0  (3.  la) 

Equation  (3.  1)  is  the  general  form  of  Bessel's  Differential  Equation  of 

order  n#  A  solution  of  (3.  1)  is 

yU)  =  Jn(x)  (3.  2) 

where  J  (x)  is  known  as  the  Bessel  function  of  the  first  kind, 
n 

We  see  that  before  we  can  derive  any  benefit  from  the  theorems  of 
section  2,  we  must  reduce  (3.  1)  to  its  normal  form.  Write  (3.  1)  in 
the  form 

i  .2 

(3.  3) 


y'l  +  1  y'  +  (1  .  £-  )  y  -  0. 

x 

X 


Section  2.  2  tells  us  that  (3.  3)  may  be  reduced  to  normal  form  by 
the  transformation 

y  =  v  exp  [  -  j  j *1  dx  ] 

=  v  exp  [  -  ^  loge  x  ] 

(3.4) 


=  v  x 


-1/2 


Then 


,  1  ,  -.3/2  ,  ,  -1/2 

y  =  -  7  v(x)  +  v1  x 


(3.5) 


y"  =  x  v  X 


-5/2  1  ,  -3/2  1  ,  -3/2  „  -1/2 

:  -  j  v'x  -  —  v’x  +  v"x 


(3.6) 


Substituting  (3.4),  (3.5),  and  (3.  6)  into  (3.3)  and  simplifying  we 


have  the  desired  normal  form  of  Bessel's  Equation 

l-4n^ 

v"  +  ( 1  +  — )  v  =  0. 

4x 


(3.  7) 


The  transformation  (3.  4)  implies  that  (3.  7)  has  a  solution 

v(x)  =  x*  ^  JJx).  (3.  8) 


n 


Trivially  any  ze ro  of  (3.2)  is  also  a  zero  of  (3.8). 
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Consider  the  differential  equation 

wM  +  w  =  0  (3.  9) 

w(a)  =  0 

We  know  that  a  solution  to  (3.  9) 
w(x)  =  k  sin(x-a) 

vanishes  at  x  =  a  and  again  at  x  =  a  +  tt  . 

Suppose  we  examine  (3.7)  and  (3.9)  with  respect  to  Theorem  2.  1 
(Sturm*  s  Comparison  Theorem).  Let  us  denote  the  coefficient  of  v  in 
(3.  7)  as 

P(x)  =  1  +  (3.  10) 

4x 

We  see  from  (3.  10)  that  we  must  consider  two  cases,  1)  when  n  =  0 
and  2)  when  n  ]>  1  . 

Case  I. 

n  =  0.  We  see  that  P(x)  >  1  and  continuous,  therefore,  theorem 
2.  1  applies.  We  conclude  that  if  the  solution  (3.  8)  to  (3.  7)  vanishes 
when  x  =  a,  it  will  vanish  again  on  the  interval  (a,  a  +  ir).  This  implies 
that  the  difference  between  consecutive  roots  of  JQ(X)  is  less  than  tt  . 
Moreover,  we  see  that 
lim  P(x)  =  1  * 

X  0=5 

This  tells  us  that  with  respect  to  the  index  i  of  roots  r^  ; 

lim  j  ri+1  -  i\  |  =  it  „  (3.  11) 

o  -> 

Case  II. 

1  >/  2 

n  ;>  1.  In  this  case  0  <  P(x)  <  1  when  /  4n  -1  <x<  &&  . 
Then  x  must  be  in  the  interval  (  ^  V  4n^  -1  ,  o*?  )  in  order  that  (3.  7) 
have  an  oscillating  solution.  Theorem  2.  1  tells  us  that  a  solution  to 
(3.  7)  which  vanishes  at  x  =  a  will  not  vanish  again  in  (a,  a  +  ir  ). 
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Therefore,  the  difference  between  the  consecutive  roots  of  J  (x)  is 

n'  8 

greater  than  *rr  .  Similarly, 

lim  P(x)  -  1,  when  n  ^  1  , 

tells  us  that  the  limiting  difference  is  the  same  as  in  Case  I,  namely  the 
difference  given  in  (3.  11)0 

The  treatment  of  the  differences  between  consecutive  roots  of  the 
Bessel  functions  in  the  above  two  cases,  suggests  that  further  treatment 
of  the  problem  also  be  made  in  two  parts,  i,  e0  1)  n  =  0  and  2)  n  J,  1  g 
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4.  Application  to  Bessel's  Equation  of  Order  Zero, 

When  n  =  0,  the  normal  form  of  Bessel's  differential  equation 
reduces  to 

v"  +  (1  +  — K-  )  V  =  0o  (4.  1) 

4x 

Suppose  the  conditions 

v(a)  =  v(b)  =  0  (4.  2) 

are  imposed.  Let  the  coefficient  of  v  in  (4.  1)  be  denoted  by  P(x)»  i.  e„ 

P(x)  =  1  +  (4.  3) 

4x 

Clearly  P(x)  y  1  for  all  x,  P(x)  is  continuous  and  decreasing,,  Then, 

from  theorem  2„  2,  there  exists  an  x  in  the  interval  (a,b)  such  that 

the  equation 

u"  +  (1  +  — L  )  u  =  0 
4  x 

has  a  solution 

u  (x)  =  k  sin  (x~a) 

with  a  root  at  x  -  a  and,  moreover,  will  have  the  next  consecutive 
root  at  x  =  b# 

Let  us  now  attempt  to  construct  a  sequence  |x^  which  will  converge 
to  x  0  There  are  two  immediate  choices  for  the  first  element  of  the 
sequence,  a  or  a  +  tt  .  We  will  generate  the  sequence  by  means 

of  the  following  recursion  formulae: 


p(x )  =  i  +  — L- 

(4.4) 

4x. 

i 

u!'  f  P(x.)  u.  =  0 
i  '  i  i 

(4.5) 

u^x)  =  k.  sin  V^jx^iix-a) 

(4.  6) 

Xi+1  =  a  +  w/Vptej) 

(4.  7) 
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Wherein  (4.  6)  is  a  solution  to  (4.  5)  which  is  forced  to  have  a  root  at 
x  =  a;  and  x^+^  is  the  next  consecutive  root  of  (4.6).  We  see  from 

(4.4)  that  P(x.)  y  1  for  all  x.  .  This  implies  that  x^+^  <  a  +  ir  for 

all  x^  . 

Case  I. 


Choose  xq  =  a,  giving 


P(xo)  -  1  +  — T  =  1  +  “7 

4x  4a 

o 

(4.8) 

u"  +  P(x  )  u  =0 
o  o  o 

(4.9) 

uQ(x)  =  kQ  sin yfP( xq) ( x- a ) 

(4.  10) 

Xj  =  a  +  it  /yf P(xQj 

(4.  11) 

Since  P(a)  P(x)  for  all  x  in  the  interval 

(a,b),  theorem  (2.1) 

(Sturm's  Comparison  Theorem)  tells  us  that 

Xj  lies  in  the  interval 

a  <  Xj  <  b  . 

Then  is  substituted  in  the  recursion  equations  giving 

p<*i> = * + ~i 

4xi 

(4.  12) 

u'j1  +  P(x^)  u^  =  0 

(4.  13) 

u^(x)  -  kj  sin  Vp(xj)  (x-a) 

(4.  14) 

x 2  —  a  +  it  /  Vp(xj) 

(4.  15) 

We  cannot  apply  theorem  2.  1  with  respect  to  (4.  13)  and  (4.3) 
since  we  cannot  guarantee  the  necessary  inequality  between  P(Xj)  and 
P(x).  However,  we  may  apply  theorem  2.  1  with  respect  to  (4.  13) 
and  (4.9)  since  P(x^)  <  P(xq),  This  tells  us  that  x^  >  Xj  . 
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If  we  expand  (40  15) 


x2  =  a 


we  see  the  sequence  is  generating  a  rather  strange  appearing  continued 
fraction;  strange  in  the  sense  that  the  more  common  continued  fractions 
do  not  contain  radical  expressions*  Moreover,  note  that  the  fractional 
form  becomes  exceedingly  more  complicated  as  each  new  term  is  added* 
Let  us  now  assume  that  for  some  m  >  0,  we  have 

xm-l  >  *m-2  (4-16> 

The  inequality  (4,  16)  implies  that  the  following  inequality  also  holds 

(4.  17) 


P(x  i)  =  1  +  -V-  <  P<x  1  +  — ~T 

1 JV  i  i  A.  ^ 

m- 1  2 


Theorem  2e  1  immediately  tells  us  that  x  is  greater  than  x 

1  m  m- 1 


Inductively  we  have  shown  that  the  sequence  |x.(  is  increasing* 
Moreover,  we  have  shown  that  the  sequence  is  bounded  above  by  a  +  it, 
consequently,  the  sequence  converges  to  some  x  such  that 

T T 


x  ~  a  + 


1  + 


1 

4x 


(v'J  \Z  ,  r "  Z  2 

x-a)  ~  4x  TT 


4x  +  1 


(4x  ^  +  1  )(x  -  a)^ 


2  2  « 
4x  v  =  0 


4x  ^  -  8a  x  %  (4a^  +  1  -  4ir^)  x  ^-2a  x  +  aC~  0  * 

Then  x  is  the  root  of  (4*  24)  that  is  nearest  a  +  it 

2 


(4„  23) 

(40  24) 

We  see  from 


(4*  23)  that  for  large  roots,  4  x  )/  1  and  the  large  roots  of  (40  24)  may 
be  approximated  by: 
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4s\*"  *4  ,  &s  .2  c\?  2  2  _ 

x  (x-a)-4x  it  =  0 

.2  2. 

(x  -  a)  -  tt  -  0 

<r*s/  o 

X  =  a  +  TT 

A  simple  calculation  will  show  for  any  given  root  !a*  of  JQ(x)  that 
(4.  24)  has  a  root  near  a  +  tt  „ 

Case  II* 

Let  us  choose  x  =  a  +  tt  * 
o 


P(x  )  =  1  + 
o 


1 


4(a  +  t t)c 


uM  +  P(x  )  u  =0 
o  o  o 


u  (x)  =  k  sin  \f  P(x  J  (x-a) 


(4.25) 

(4.26) 


x^  =  a  +  tt  / 


xo) 


(4.2  7) 

(4.  2  8) 

We  have  P(a  +  it)  <  P(x)  for  all  x  in  (a,b),  therefore,  we  know  x 
cannot  he  in  (a,b),  therefore,  we  see 


1 


b  <  Xj  <  a  +  ir 

Substituting  in  the  recursion  formulae  we  have 

1 


P(Xj)  =  1  + 


4x 


1 


Uj"  +  P(xj)  nl 


u^x)  =  kj  sin  yfP(x1)  (x-a) 
x.^  =  a  +  t r  /  VPU.) 


(4.  29) 

(4.  30) 

(4.  31) 
(4.32) 
(4.  33) 


The  inequality  ^4.29)  tells  us  that  P(x^)  y  P(xc),  and  we  get  x^  <  x^ 
from  Theorem  2.  1.  Moreover,  (4.  29)  also  tells  us  that  P(x^)  <  P(x) 

for  all  x  in  (a,b)  and  x^  b  . 

If  we  reverse  the  assumed  inequality  in  (4.  16)  we  can  show  induc¬ 
tively  that  this  sequence  £x.f  ,  generated  by  starting  at  xq  -  a  +  tt  , 
is  a  monotone  decreasing  sequence.  Moreover,  we  see  that  ,|x^C  >  b 
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for  all  i,  hence  is  bounded  below*  Then  the  sequence  converges  to  an 
x  j,  where 

x j  =  a  +  it  /  s[p('kl  ) 


TT 


(4  +  l)(xj-a)*  *  4  tt  =  0 

4  xx4  -  8a  +  (4a^  +  1  -  4tr  )  x^  -  2a  -  0  (4*  34)  0 

The  quartic  (4*  34)  has  coefficients  which  are  identical  to  the  coef¬ 
ficients  of  like  powers  of  x  in  (4*24)  ,  therefore  we  conclude  that 
Xj  =  x  *  This  shows  that  there  is  no  difference  in  the  limit  between 
the  sequence  generated  in  Case  I  and  that  generated  in  Case  II*  Un- 
fortunately,  we  have  also  shown  that  x  y  b  ,  but  we  have  no  indication 
of  the  difference  between  x  and  b* 

The  quartic  (4*  24)  has  been  solved  with  the  aid  of  the  CDC  1604 
Digital  Computer  for  the  first  20  roots  of  Jo(x)0  The  computed  values 
of  x  are  compared  with  the  tabulated  values  of  the  roots  of  J^(x)  «  The 
difference  between  x  and  the  tabulated  value  of  b  is  shown  in  the  error 

column*  The  difference  has  its  maximum  at  the  first  root  and  decreases 
quite  rapidly  as  x  grows  larger* 

Simple  calculations  show  that  xT  is  almost  midway  between  b  and 

a  +  tt  for  the  first  root,  but  this  relation  does  not  persist  as  the  root 
number  increases,  In  fact  x  approaches  b  more  rapidly  than  a  +  tt  „ 

The  fact  that  the  sequence  x^  does  not  converge  to  b  implies  that 

some  change  to  the  recursion  formulae  must  be  made,  The  most  direct 

method  of  changing  the  formulae  consists  of  a  modification  to  P(x)w 
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4.  lv  Linear  Average* 


Suppose  a  sequence 


were  to  converge  such  that 


(4.  1.  1)  holds t 


i.  e. 


b 


a  +  tr 


where  & 


is  the  limit  point  of  a  sequence 
(4%  2  +  l)(b-a)2  =  4%  2  tt2 

4%  2  [  (b-a)2  -  tt2  ]  =  -  (b-a) 


4 


(4.  1.  1) 


We  have  shown  that  (b-a)  <  tt  ,  therefore,  we  know  that  %  is  real* 

Computed  values  of  (  are  displayed  in  table  II ,  along  with  the  tabulated 

values  of  b.  The  error  column  is  the  difference  between  x  and  %  , 

and  indicates  the  amount  x  must  be  decreased  if  the  sequence  x/  is  to 
converge  to  b.  Noting  that  %  is  almost  midway  between  a  and  x  the 

following  modification  to  the  recursion  formulae  is  suggested: 


1 


P(x.) 

=  1  +  2 

4(a  +  x.) 

(4.  1.2) 

— Z - 

u.  + 

1 

P(x.)  u*  =  0 
i  i 

(4.  1.3) 

vu(x) 

=  sinVp(x^)  (x-a) 

(4.  1.4) 

xi+l 

=  a  4  tt  /  /p(x.) 

(4.  1.5) 

In  essence,  the  modification  consists  of  merely  substituting 
(a  4  x  ) 

'  r  in  place  of  x.  in  (404)«  Clearly  the  conclusions  reached 

— F“  1 

<s-  y 

regarding  convergence  of  the  sequence  ; ’  are  also  valid  for  the 
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sequence  f  „  For  the  sequence  y  we  are  unable  to  make  any  con¬ 
clusion  relative  to  b  since  the  of  Ca.se  II  is 


x 


1 


a  +  tt  / 


a 


+ 


i r 


"7 

) 


(4.  1.6) 


and  lies  in  the  interval  (a,b)e  Moreover,  we  can  see  that  for  i  ^  1 
each  element  of  the  sequence  f  x^|  is  smaller  than  the  corresponding 


element  of  the  sequence  and  we  conclude  that  the  limit  of  the  se¬ 


quence 


fxit  * 


x  is  smaller  than  x  .  We  have  that 


x  =  a  + 


'\AT=I 


=  a  + 


,/a+x» 

M-T-) 


VI  + 


(a  +  x) 


(x  -  a)2[  (a  +  x)2  +  I  ]  -  (a  +  x)2  tt2  =  0  (4.  1.  7) 

x  +  (1  -  2a2  -  x2)  x  2-  [  2a(  1  +  tt2)]  x  +  a4  +  a2(l-ir2)  =  0  (4.  1.  8) 

For  the  larger  roots  of  (4W  lQ  8)  we  look  at  (4*  I„  7)  and  see  that 
_  2 

(a  +  x)  y  y  1  so  that  an  approximation  yields 
(x  -  a)"  (a  +  x)^  -  (a  i  x)^  ir^  =  0 
x  =  a  +  it  0 

We  conclude  that  (4*  1„  8)  has  a  solution  near  a  +  tt,  and  that  this 
solution  is  the  x  we  desire0 

Table  III  lists  the  desired  roots  of  (4.  1*8)  for  the  same  20  roots 
as  in  section  4W  A  comparison  of  the  error  column  of  table  III  against 
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table  I  shows  that  the  difference  between  x  and  the  tabulated  root  b  is 


about  one-tenth  as  large  as  the  difference  between  x  and  b  at  the  first 


root.  Moreover,  we  see  that  x  approaches  b  more  rapidly  than  x, 
that  x  -  b  is  only  .  01  times  K  -  b  for  the  tenth  root. 

B  ased  upon  the  results  of  these  computations,  we  conclude  that 
although  x  is  nearer  b  than  x  ,  x  is  still  not  the  value  of  x  whose 
existence  is  guaranteed  by  theorem  2.  2. 

4.2.  Mean  Value. 

Consider  the  integrated  mean  value  of  P(x)  over  the  interval  (a,x) 

and  denote  this  integrated  mean  value  by  F(x),  then 

t  x  , 


F(x)  = 


=  —  f 

x"a  J  a. 


(1  + 


4t 


)  dt 


1 


x-a 

1 _ 

x-a 


Ft  -  —  lx 

1  4t  Ja 

r  ,  1  /  1  1m 

[x-a  +  t  (  -  -  -  )  J 


1 


=  1+  4^r  ) 


=  i  + 


4ax 


We  see  that 


1P(x)  =  1  +  4^  >  P(X)  =  1  +  — ~ 

4ax  (a+x)^ 

i  s  1 

4ax  '  (a+x)2 

a2  +  2ax  +  x2  ^  4ax 
a2  -  2ax  +  x2>  0 
(a-x)2  >  0  for  all  x  4 

Suppose  we  generate  a  sequence  x.‘ 

X  <> 

recursion  formulae: 


(4.2.1) 


a  . 

by  means  of  the  following 
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lP(x.)  =  1  +  - 

i  4ax. 

x 

(4.  2.2) 

u.M  +  IP(x. )  u.  =  0 
i  1'  l 

(4.2.3) 

u.(x)  =  k.  sin  "V  lP(x.’)(x-a) 

(4.2.4) 

Xi+1  =  a  +  wVlPU.) 

(4.  2.5) 

But  equations  (40  2S  2)  through  (4.  2.  5)  are  the  same  as  the  set  of 

2 

equations  (4.4)  through  (4.7)  except  ax.  is  substituted  for  in 

(4.4).  If  we  make  this  substitution  throughout  our  preceeding  considera- 

•  C  1  A 

tion  we  see  that  the  sequence  |  x.  also  converges  to  some  x  =  x  such 

that 


(4ax  +  l)(x-a)^  -  4ax  tt^  =  0 

4ax 3+  (l-8aZ)  x2  +  (4a3  -  2a  -  4a  it2)  x  +  a2  =  0  (4.  2.9) 

Inequality  (4.  2.  1)  shows  us  that  x  must  indeed  be  less  than  x 
but  once  more  we  are  unable  to  state  any  analytical  conclusions  regard¬ 
ing  x  and  b  . 

Table  IV  contains  the  roots  of  (4.  2.  9)  for  the  same  zeros  of 
JQ(x)  as  previously  computed.  We  see  that  x  is  smaller  than  b.  We 
note  also  that  |x-bj  is  greater  than  |x-b|  . 

Since  b  has  been  bracketed  by  x  and  x  *  a  linear  average  of  the 
two  is  suggested.  Accordingly  we  define 

P(xl.)  a  (P(x.)  +  P(x.))/2  (4.2.10) 

xl.  +  i  =  a  +  it/  Yp(x1.)  (4.  2.  11) 

Then  equation  (4.  2.  10)  and  (4.  2.  11)  together  with  equations  (4.  5)  and 
(4.  6)  constitute  a  set  of  recursion  formulae  which  will  generate  a 
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sequence  fxl  .  We  know  the  sequence  xl  will  converge  to  some 

_  r  t  ' 

xl  since  it  is  the  sum  of  two  convergent  sequences  ami  x^  .  The 

values  of  xl  are  displayed  in  Table  V.  We  see  that  although  the  values 
of  xl  are  nearer  to  b  than  x,  they  are  not  as  near  as  x  . 

Suppose  we  consider  the  average  of  x  and  x  .  Let  us  define: 

P(xZ.)  =  (P(x.)  +  P(x.))/2  (4.2.12) 

x2i+l  =  a  +  VP(x2.  )  (4.2.13) 

Again  equations  (4.2.12)  and  (4.2.13)  together  with  (4.5)  and  (4.6) 

will  generate  the  convergent  sequence  x2.  .  We  see  from  Table  VI 

that  x2  is  nearer  to  b  than  any  of  the  previous  values  generated*  The 

-3 

maximum  difference  between  x2  and  b  is  less  than  10  for  the  first 

~  7 

root  and  decreases  to  less  than  10  for  the  20th  root.  The  difference 
decreases  quite  rapidly  as  is  evidenced  by  3.  value  less  than  10  ^  beyond 
the  fifth  root. 

4.3  Cumulative  Effect. 

Let  us  now  investigate  the  cumulative  error  for  each  of  the  sequen¬ 
ces  previously  discussed.  The  errors  in  this  case  are  cumulative  in 
the  following  sense: 

Assume  that  a  =  r^  the  first  root  of  JQ(X)  is  known,  we  will 
compute  an  x  from  a  sequence.  Then  each  successive  root  is  generated 
using  the  root  last  computed  as  the  value  of  the  known  root. 

The  cumulative  errors  for  x,  x  ,  x,  xl  and  x2  are  shown  in 
Tables  VII  through  Xi  respectively.  As  well  expected,  the  cumulative 
error  is  due  mostly  to  the  error  in  computing  the  first  root.  We  see  from 
Table  XI  tha,t  the  cumulative  error  for  the  20th  root  is  only  10% 
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greater  than  the  error  for  the  first  root.  The  worst  case  is  shown  in 
Table  X  where  the  cummulative  error  for  the  20th  root  is  40%  greater 
than  the  error  for  the  first  root* 
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(5.  1) 


50  Application  to  Bessel1  s  Equation  of  Order  n  ;>  1  . 

When  n  J>.  1,  the  normal  form  of  Bessel's  equation  is 

y"  +  U  +  )  y  =  o 

4x 

Let  us  impose  the  condition 
y(a)  =  y(b)  =  0  . 

Again  let  P(x)  denote  the  coefficient  of  y  in  (5.  1),  i.e. 

1  -4n^ 

P(x)  =  1  +  (5.  2) 

4x 

We  see  that  P(x)  >  0  for  all  x  y  ^  v^rr,  moreover,  P(x)  is 
monotone  increasing,  lim  P(x)  =  1.  i.  e0  0  <  P(x)  <  1  for  all  x 

x  ^  i  j — pT — 

such  that  x  is  in  the  interval  ^  y4n  -1,  ^  )  .  We  shall  limit  our  dis¬ 
cussion  to  the  case  P(x)  >  0  in  order  that  the  conditions  of  the  theorems 
of  section  2  may  be  satisfied* 

Except  for  P(x.),  the  recursion  formulae  defining  the  sequence  will 
be  the  same  as  in  section  4,  i„e. 


P(x  )  =  1  + 

(5.  3) 

4x 

o 

u!’  +  P(x.)  u.  =  0 
i  i  i 

(5.4) 

u^(x)  =  k;  sinVP(x^)  (x-a) 

(5.5) 

x.  .  =  a  +  Vp(x.) 

l+l  y  i 

(5.6) 

As  before  we  consider  the  two  cases  x  =  a  and  x  =  a  +  tt  . 

o  o 

Case  I: 

Choose  x  =  a  ,  giving 

o  o  o 


tp i/  \  i  ,  l-4n  _  ,  l-4n 

p(xo)  =  1  +  „  TT  "  1  +  TT 

4x  4a 

o 

(5.  7) 

uM  +  P(x  )  u  =0 
o  o  o 

(5.8) 

uq(x)  =•  kQ  sin  V  P(xq)  (x-a) 

(5.9) 

x1  =  a  +  it  /  VP(x0) 

(5.  10) 
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Since  P(x)  is  increasing,  we  see  that  P(xq)  <  P(x)  for  all  x  in 

(a,b)w  Then  theorem  2.  1  tells  us  that  b  .  Moreover,  P(xq)  <  1  > 

therefore,  from  (5*  10)  we  see  that  >  a  +  tt  . 

When  x,  is  substituted  in  the  recursion  equation  we  have 
1  2 

P(x  )  =  1+  l-^p-  (5.11) 

4xx 

Uj  +  P(x^)  u^  =  0 
u^(x)  =  sin  Vp(  )  (x-a) 

x-,  =  a  +  tr  /^P(xj)  (5.  1Z) 

We  know  from  the  inequality  x^^b  that  x^  is  in  the  interval  (a,b). 
Trivially  then  x^  <  x^,  moreover,  P(x^)  <  1  tells  us  that  x^  >  a  +  *rr  . 
Then  we  have 


x  =a<a+7r<x:)<b<x1 
o  2  1 


(5.  13) 


Let  us  now  consider  the  next  element  of  the  sequence.  If  x^  is 


substituted  in  equations  (50  2)  through  (5.  6)  ,we  have 
P(x  )  =  1  +  h|*i 

4xz 


(5.  13) 


u 


2  +  P(xZ)  u2  =  ° 


u,(x)  -  k?  sin  V  P(x?)  (x-a) 


x^  =  a  +  tt/Vp(x2) 


‘2'~'  "  ~2  '  2' 

(5.  14) 

Since  is  in  the  interval  (a,b)  we  do  not  have  the  necessary 

inequality  for  theorem  2*  1  and  consequently  we  can  no  longer  determine 
the  relation  of  x^  to  b  for  i  J>  3  „  We  see,  however,  that  P(x^  )>P(Xj) 
since  P(x)  is  increasing  and  x^  <  .  This  implies  that  x3  ^  x2  • 

Clearly  P(x^)  <  P(xq)  and  x^  <  „  Then  we  have 


x  s  a  <  a  +  u  <  x-,  <  x0  <  x1 
o  c,  />  l 


(5.  15) 
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Suppose,  for  some  k^3  that  the  following  inequalities  hold: 


a  <  a  f  Tr<Xm  <  «  •  •  <  x,  <  x,  <  x,  < 
^  k-  2  k  k  **  I 


<  x 


The  recursion  formula  gives 

p(*k  ?)  *  1+  i“4n" 


k-2 
2 


p(xk) 


l  + 


P(xk_  j)  =  1  + 


4  x 

1-4n 

-T" 

4ltk 

i  „  i 

l  -4n 


(5.  16)  . 

(5.  17) 

(5.  18) 

(5.  19) 


4x. 


■k-1 


uk-2  +  P(xk-2)uk-2  s  0 

(5,  20) 

+  P(xk)uk  -  0 

(5.21) 

uk- 1  +  P*xk-Puk-i  “  0 

(5.22) 

xk  1  =  a  4- -it  /  Vp(x_k  p)  (x-a) 

(5.23) 

xk  =  a  +  tr  /  v  P(xk  t )  (x-a) 

(5.  24) 

xk+l  "  a  +  17  ^  P(xk)  (x-a) 

(5.  25) 

The  assumption  (5.  16)  tells  us  that  P{x^  >)  ^P(xk)^>P(xk  )  since 
P(x)  is  increasing.  If  we  apply  theorem  2.  ls  first  with  respect  to  (5.  21) 
and  (5,20),  then  with  respect  to  (5„  2'i)  and  (5.22);  we  have 

xq  =  a  <  a  f  tr  <  <  „  . .  <  xk  ^  <  xk  <  xkk^<  xk  ^ <„  .  „  <x^  (5.  26) 

We  have  now  shown  inductively  that  ..  lies  in  the  interval 


(x. ,  x.  ,)  for  all  i  .  Let  us  now  select  the  sub  sequence  -fz  A  of  |x.  I 

1  i-  I  ^  *  1 9 

such  that  for  each  j;  z  .  ~  0  Trivially  z  i  25  ^  z  =  x  0  Assuming 

j  r  o  o 

for  some  m  >  0  wf  see  that  z 


then  z  >  z 
m 


m 


x2m>x2m-2IIm-r 


A  m- 1 

Moreover^  we  know  that  x  ^  lies  in  the  interval  (x-,  ,  %  .  ,  ) 

v  •  m  i  i  tj  m  &*>  m  t  I 


therefore,  we  conclude  z 


mil 


y  z 


m 


and  that  fz  V’  is  monotone  increasing. 
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Clearly*  the  sequence  z..  ^xs  bounded  abov?s  by  *  «  Therefore*  we 
conclude  that  the  sequence  z.  converges  ro  some  z  „ 


Similarly*  consider  the  subsequence  |vn|  of  -'.'x.  chosen  such  that 

T'  7/ 

for  each  n*  s  1°  an  argument  similar  to  the  one  for  z ^  we 

can  show  that  the  subsequence  is  decreasing  and  bounded  below  by 

a*  therefore*  it  converges  to  some  v  0  We  have  now  shown  that  the 
sequence  | x/|  contains  two  convergent  subsequences  and  v^  *  and 
in  order  to  prove  convergence  of  |  we  must  show  that  the  subsequences 
converge  to  the  same  point*  i0  e0  that  z  ^  v0  We  know  from  (5„25) 
that 


a  4' 


IT 


a.  f 


TT 

(7 

l-4n2 

4  /  a  + 

(5.  26) 


ir 

aTeS! 

f  4  7  ‘ 


and  v 


a  + 


/ 


l-4n 

4  /  •  t  ,.4= 


\ 


1  t 


l-4n 

7^“2  i 
4v  / 


(b.  2  7) 


Equation  (5.  26),  after  considerable  algebraic  manipulation,  can  be 
reduced  to  the  equation 

[z-a)^(4  l  a(4z  f  i-4n‘  )  f  4z  +  (l-4n  }(4z  +  l-4n‘ ' )-» 4-ir" 

£a(4z  -f  I“4n  )  ^  4z  tr  |  ]  I6a^5s~(4z;  ir  ,  (z^  a)”-  4?r  ^  0o  (t>„  2{ 

which  is  clearly  a  polynomial  of  degree  8  M  Similarly*  Equation  (5„  27) 
will  reduce  to  the  same  polynomial,  Then  we  can  conclude  that  z  -  v 


and  the  sequence  1>-  | converges  to  some  x  such  that 
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x  =  a  + 


which  reduces  to 


IT 


1  + 


1-  4n 


4x 


4x  ^-8ax  ^  +  (4a  "+  l-4n  ‘  ~4ir^)x  2a(  l-4n^)  x 

+  (l-4n2)  a2  =  0  (5.  29)  . 

Moreover,  x  is  the  root  of  equation  (5C  29)  nearest  a  +  tt  . 

Case  II: 

Suppose  we  choose  x  =  a  +  tt  then  (50  3)  through  (5.  6)  becomes 

P(x  )  =  1  +  Ills  =  1  +  ills.  ?  (5.  28) 

4x.  4(a+Tr) 

u"  +  P(x  )  u  =  0 
o  o  o 


uq(x)  =  kQ  sinV^P(xQ)  (x-a) 
x^  =  a  +  tt  /  Vp(^0) 


(5.29) 


We  have  chosen  an  xq  with  interval  (a,b)  ,  therefore,  we  cannot 
make  any  comparisons  between  x^  and  be  Since 
P(x  )  <  1  we  know  x,  >  a  +  tt  =  x 
Let  us  now  find  the  next  element  x^  of  the  sequence  ^  x.|  .  We 
have  the  equations 


P(Xi)  =  j  +  Id* 


2 


4x 


1 


u'j*  +  P(.Vj)  Uj  =  0 
Uj(x)  -  k,  sin  Vp(x j)  (x-a.) 


The  inequality  x.  >  >  implies  that  P(x  J  <  P(x  )  since  P(x)  is 


9  therefore,  tells  us  that  , 


decreasing.  Theorem  2, 
moreover,  from  P(> )  <  1  we  see  that  x^  ^  a.  +  tt  giving 


x 


a  +  it  <  X£  <  Xj 


(5.  30) 
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Clearly,  ior  all  i>  2  this  sequence  behaves  exactly  as  the  sequence 


ol  Case  I,  therefore,  converges  to  somt'  Xj  such 


that 


x^  -  a  + 


ir 


1  + 


i  -4n 


4^ 

X 


(5.31) 

But  equation  (503Q1)  after  simplification  reduces  to  the  same  quartic  as 

.  vV  if'i 0 

(5«  2  7),  and  we  conclude  that  =  x  0 

In  either  case  for  n  >1  we  have  been  unable  to  relate  x  to  b, 
therefore,  we  must  consider  computations  for  the  first  20  roots. 

The  results  of  the  computation  displayed  in  table  XII  are  com¬ 
parable  to  those  shown  in  table  T  for  the  case  n  s  0,  as  should  well  be 

expected.  We  see  that  for  the  20th  root,  the  error  in  table  XII  is 

-4  -  6 

10  while  for  J  ,  it  is  about  5  x  10 

o 

Suppose,  however,  the  sequence  w *•  r*  to  converge  to  b.  The 

final  result  would  for  some  %  ,  the  following 

IT 

b  =  a  f 

'PU  ) 


=  a  + 


1  + 


1  -  4n 


(4X  2  i  1  4 r,7)  {b- a)“  -  4  i  2  tt 

which  when  solvec  (or  %  snows  that 

K&  ^  ( ’5-“  a ) 


((b-a)2-  it2) 


(5.  32) 


Trivially  f  determined  in  this  manner  is  real  valued,  since 
(b-a)  ^  it  and  4n2  y  1  „  The  values  of  %  computed  from  equation  (5.32) 
are  shown  in  table  XIII.  These  values  suggest  fha?  the  following  modifi¬ 
cation  be  made  to  the  recursion  formulae. 
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PC*,) 


i  -fr 


j__  4n 

„  f v  . . 4 
i) 


u  "  +  P(x  )u  =  0 

1  3  3 


u.(x)  s  k.  sin  fP(x.)(x-a) 


l  i  '  r 

x. .  .  =  a  + 


IT 


i+1 


VP(x.) 


We  have  already  shown  that  the  sequence  l)x.  defined  by  equations 


(5.3),  (5.4)s  (5.5)  and  J5„  6)  converges  to  x.„  The  sequence  •  '  is 

the  sum  of  a  constant  and  a  multiple  of  x..,  and  therefore  will  converge. 
This  implies  that  the  sequence  |" converges  to  some  x  such  that 


x  =  a  + 


/I  t 


Un_ 

M^) 


(5.  13) 

y  4,;~r) 

x^+(l-4n^-2a^-ir  ')x*  -  [ 2a(J ~4n^4ir^)  ]x  4  a4f  a  ~ir~)  -  0  (5.  14) 

The  roots  of  the  quar&c  (50  14)  nearest  a  r  nr  are  displayed  in 
table  XIV  for  the  first  20  roots  of  J^(x)w  Comparing  the  errors  of 
ta-ble  XIV  against  table  XII  we  see  that  Jx~bj  is  one  order  of  magni¬ 
tude  smaller  than  ]  *  b  j  „ 

We  shall  nex<  .-amine  the  sequence  |  \  ,  i0e0  the  mean  valuer  of 

P(x)  over  the  int<  rv.*l  (a,  The  mean  value  o’  P{\)  is  defined  in  the 
same  manner  as  in  section  4„  20  We  know  that  the  sequence  x.  f  will 
converge  since  it  consists  of  elements  which  are  multiples  of  the  elements 
of  a  convergent  sequence tt  Then  we  have 


x  =  a  +  it/  ¥  P{y 

ir 


a  4* 


i 


i  + 


•  4n 
4  a  x 
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(5.  15) 


4axV(l-4n^- 8a  )x^+{4a^-2a*8a^ni -4anr)x  a  -  0 

We  see  that  for  n  £  1,  the  limit  of  the  sequence  results  in  a  cubic 
similar  to  the  one  for  x  when  n  =  0*  Again  we  know  that  x  is  the  root 
of  the  cubic  nearest  a  +  tt. 

Table  XV  displays  the  roots  of  equation  (5*  15)  for  the  first  20 
roots  of  Jj(x)„  As  expected  from  the  results  for  JQ(X)»  we  see  that 
again  x  is  past  b  and  that  x  and  x  have  straddled  b.  We  note, 
as  for  Jq(x),  }x>bj  >  |x-b|  . 

We  shall  now  examine  xl  and  x2,  as  defined  in  section  4.  2.  The 

values  of  xl  are  displayed  in  table  XVI  „  Again*  as  for  Jq(x),  we  see 

that  xl  is  further  from  b  than  x„ 

x2  is  again  nearest  b  for  J^(x),  with  a  maximum  error  for  the 

first  root.  In  general,  we  see  that  the  various  sequences  generated  for 

Jj(x)  tend  to  converge  to  points  farther  from  D  than  the  same  sequence 

for  J  (x). 
o 

As  expected,  the  cumulative  errors  shown  in  tables  XVIII  through 
XXII  exhibit  the  same  characteristics  as  the  cumulative  errors  for 
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6.  Conclusions. 


We  have  proposed  a  technique  of  computing  the  consecutive  zeros 
of  a  solution  to  a  Sturm- Liouville  System. 

[r(x)y'(x)]'  +  P(x)y(x)  =  0  (6.1) 

y(a)  =  0  y(b)  =  0 

The  technique  requires  the  determination  of  an  x  such  that  the  equation 

[r(x)u'(x)]'  +  P(x)y(x.)  =  0  (6.2) 

will  have  a  solution  u(x)  which  maybe  forced  to  vanish  at  x  =  a  and 

at  x  =  b  . 

The  existence  of  such  an  x  has  been  shown  by  theorem  2.2.  More¬ 
over,  in  theorem  2.  2.  1  we  have  established  upper  and  lower  bounds  on 

«*v 

the  range  of  this  x  . 

In  the  specific  case  of  the  Bessel  functions,  we  note  that  each  of  the 

sequences  generated  by  the  proposed  technique  converges  to  a  particular 

x,  no  matter  whether  the  first  element  was  chosen  equal  to  the  known 

root  a,  or  chosen  equal  to  a  +  it  .  We  see  that  x  turns  out  to  be  the 

root  of  a  quartic  nearest  a  +  it.  Similarly,, we  see  that  x  is  also  a  root 

A 

of  a  quartic  polynomial  while  x  resolves  into  a  cubic  equation. 

The  solution  for  %  for  the  case  n  =  0  and  for  n  =  1  both  indicate 
that  the  value  of  t\v  desired  x  as  defined  by  theorem  2.  2  is  almost 
midway  between  a  and  b0  The  difference  between  x  and  %  approaches 
tt/2  as  x  increases* 

We  have  also  shown  that  the  desired  x  for  the  Bessel  functions  lies 
between  the  x  determined  by  a  linear  averaging  technique  and  the  x 
determined  by  the  mean  value  of  P(x)  over  the  interval  (a,  x)0 

Unfortunately,  all  of  the  above  generating  techniques  failed  to  result 
in  a  sequence  converging  to  the  desired  x„ 
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TABLE  I 


VALUES  OF  X  ***  JqIXJ 


ZEROS 

OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO 

FOR  ROOTS 

ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION 

FORMULAE 

pixiin  =  i  +  i/Mxmr 

X(I+1>  * 

a +7*7  v>(xm ) 

ERROR 

IS  THE  DIFFERENCE 

BETWEEN  TABULATED 

AND  COMPUTED  ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

5.5200781103 

5.5336719233 

.0135938131 

2 

8.6537279126 

8. 6564432415 

.0027153292 

3 

1.1  .7915344390 

1 1.7925004759 

.0009660370 

4 

14.9309177082 

14.9313671612 

.0004494532 

5 

1 8.0710639674 

18.0713085639 

.0002445965 

6 

21.2116366294 

21.21 17842031 

.0001475737 

7 

24.3524715304 

24.3525673207 

.0000957905 

8 

27.4934791317 

27.4935447974' 

.0000656658 

9 

30.6346064676 

30.6346534281 

.0000469606 

10 

33.7758202134 

33.7758549489 

.0000347356 

11 

36.9170983527 

36.9171247659 

.0000264135 

12 

40.0584257636 

40.0584463133 

.0000205500 

13 

43.1997917127 

43. 1998080136 

.0000163012 

14 

'46.3411883712 

46.3412015187 

.0000131483 

15 

49.4826098960 

49.4826206546 

.0000107595 

16 

52.6240518400 

52.  6240607537 

.0000089143 

17 

55.7655107528 

55.7655182229 

.0000074702 

18 

58.9069839241 

58. 9069902431 

.0000063197 

19 

62.0484691886 

62.0484745828 

.0000053946 

20 

65. 1899647992 

65. 1899694391 

.0000046415 

i 
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TABLE  II 

VALUES  OF  ***  J0(X) 


ZEROS 

OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO 

FOR  ROOTS 

ONE  THROUGH  TWENTY,  BASED 

ON  THE  RECURSION 

FORMULAE 

CHI 

IS  THE  VALUE  X  MUST 

ASSUME  TO  CONVERGE  TO  B 

50  =  ( 18- 

■A  )/2  )  C  Vl/<  (B-A)2- 

-7T2)  ) 

ERROR 

IS  THE  DIFFERENCE 

BETWEEN  X  AND  CHI 

ROOT 

NO, 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

3.8368678503 

1.6968040731 

2 

8.6537279126 

7.0180586502 

1.6383845913 

3 

11.7915344390 

10. 1751194785 

1.6173809974 

4 

14.9309177082 

13.  3249362609 

1.6064309003 

5 

18.0710639674 

16.4716300704 

1.5996784936 

6 

21.2116366294 

19.6166964039 

1.5950877941 

7 

24.3524715304 

22. 7608070695 

1.5917602512 

8 

27.493479131 1 

25.9043037482 

1.5892410493 

9 

30.6346064676 

29.0473744823 

1.5872789458 

10 

33.7758202134 

32. 1902150838 

1.5856398651 

11 

36.9170983527 

35.3327245601 

1.5844002059 

12 

40.0584257636 

38.4751889762 

1.5832573371 

.  13 

43.1997917127 

41.6175294435 

1.5822785702 

14 

46.3411833712 

44.7596824104 

1.5815191083 

15 

49.4826098960 

47. 9016739726 

1.5809466820 

16 

52.6240518400 

51.0439078268 

1.5801529270 

17 

55.7655107528 

54.  1856175233 

1.5799006997 

18 

58.9069839241 

57.3278564233 

1.5791338198 

19 

62.0484691886 

60.4698329102 

1.5786416726 

20 

65.1899647992 

63.6115127150 

1.5784567241 
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TABLE  III 

VALUES  OF  X  #**  Jq(X) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 
P(X(I))  =  1+1/4(  (  (A+7(I  )  )/2)?  X(I  +  U  =  A+7T/  \/p{7{  I ) ) 


ERROR 

IS  THE  DIFFERENCE 

BETWEEN  TABULATED 

ANO  COMPUTED  ROOTS 

ROOT 

NO. 

TABULATED 

ROCT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5. 5217120022 

.0016338920 

2 

8.6537279126 

8.6538810786 

.0001531661 

3 

11.7915344390 

1 1. 7915695005 

.0000350617 

4 

14.9309177082 

14. 9309296827 

.0000119747 

5 

18.0710639674 

18. 0710691069 

.0000051396 

6 

21.2116366294 

21.2116391868 

.0000025579 

T 

24.3524715304 

24. 3524729423 

.0000014123 

8 

27.4934791317 

27.4934799736 

.0000008423 

9 

30.6346064676 

30. 6346070007 

.0000005335 

10 

33.7758202134 

33.7758205654 

. .0000003527 

1  1 

36.9170983527 

36.9170985967 

.0000002443 

12 

40.0584257636 

40.0584259368 

.0000001732 

13 

43.1997917127 

43. 1997918384 

.0000001260 

.  14 

46.3411883712 

46.3411884652 

•0000000946 

15 

49.4826098960 

49. 4826099686 

.0000000730 

16 

52.6240516400 

52.6240518950 

.0000000556 

17 

55.7655107528 

55. 7655107975 

.0000000452 

18 

58.9069839241 

58.9069839586 

.0000000349 

19 

62.0484691886 

62. 0484692166 

.0000000280 

20 

65.1899647992 

.  65.  18996482 1 6 

.0000000226 
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TABLE  IV 

VALUES  OF  $  ***  Jq { X » 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULA 
P($(I)J  =  1+1/4A$(I)  $(I  +  1)  *  A+T /VIM 

ERROR  IS  THE  DIFFERENCE  BETWEEN  TABULATED  AND  COMPUTED  ROOTS 


k 


x<  i  n 


ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5. 517232S156 

-.0028455946 

2 

8.6537279126 

8.6534819095 

-.0002460030 

3 

11.7915344390 

11.7914791510 

-.0000552879 

4 

14.9309177082 

14.9308989581 

-.0000187501 

5 

18.0710639674 

18. 0710559469 

-.0000080201 

6 

21.2116366294 

21.21 16326448 

-.0000039842 

7 

24.3524715304 

24.3524693325 

-.0000021976 

8 

27.493479 131 t 

27.4934778223 

-.0000013094 

9 

30.6346064676 

30.6346056401 

-.0000008275 

10 

33.7758202134 

33.7758196630 

-.0000005496 

11 

36.9170983527 

36.9170979746 

-.0000003774 

12 

40.0584257636 

40.0584254945 

-.0000002689 

.  13 

43.1997917127 

43.1997915152 

-.0000001969 

14 

46.341 1883712 

46.3411882240 

-.0000001468 

15 

49.4326098960 

49.4826097842 

-.0000001110 

16 

52.6240518400 

52.6240517525 

-.0000000871 

17 

55.7655107528 

55.7655106848 

-.0000000672 

18 

58.9069839241. 

58.9069838692 

-.0000000548 

19 

62.0484691886 

62.0484691439 

-.0000000444 

20 

65.1899647992 

65. 1899647620 

-.0000000366 
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TABLE  V 

VALUES  OF  in  ***  J0 ( X ) 

ZEROS  OF  THE  BESSEL  FUNCTION  Op  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE _ 

ptxfm)  =  (P(xm)+P($m))/2  TTu+i)*  a+b?  VpiTi  i  i  > ) 

ERROR  IS  THE  DIFFERENCE  BETWEEN  TABULATED  AND  COMPUTED  ROOTS 


ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5254224626 

.0053443524 

2 

8.6537279126 

8.6549613329 

.  00 12334206 

3 

11.7915344390 

11. 7919896499 

.0004552110 

4 

14.9309177082 

14.9311330232 

.0002153153 

5 

18.0710639674 

18.0711822445 

.0001182773 

6 

21.2116366294 

21.2117084200 

.0000717907 

7 

24.3524715304 

24.3525183247 

.0000467948 

8 

27.4934791317 

27.4935113089 

.0000321774 

9 

30.6346064676 

30. 6346295336 

.0000230662 

10 

33.7758202134 

33.7758373059 

.0000170928 

1  1 

36.9170983527 

36.9171113698 

.0000130179 

12 

40.0584257636 

40.0584359039 

.0000101405 

13 

43.1997917127 

43.  1997997640 

•0000080521 

14 

46.341 1883712 

46.3411948718 

.0000065007 

15 

49.4826098960 

49»  4826 152194 

.0000053242 

16 

52.6240518400 

52.6240562536 

.0000044135 

17 

55.765L107528 

55.7655144539 

.0000037015 

18 

58.9069839241 

58.9069870561 

.0000031325 

19 

62.0484691886 

62.04847186  34 

.0000026751 

20 

65.1899647992 

65. 1899671014 

.0000023024 
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TABLE  VI 


VALUES  OF  X2  ***  Jq { X ) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  OROER  ZERO  FOR  ROOTS 

ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE _ 

pfx2(in  ={pfxm  )+p(x(ij  n/2  x2(i  +  n  *  i+tyV pi xii n ) 

ERROR  IS  THE  DIFFERENCE  BETWEEN  TABULATED  AND  COMPUTEO  ROOTS 


ROOT 

NO* 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5194687883 

•".000609321  9 

2 

8.6537279126 

8. 6536814596 

-.0000464527 

3 

11.7915344390 

1  1.7915243239 

-.0000101150 

4 

14.9309177082 

14.9309143201 

-.0000033879 

5 

18.0710639674 

18.0710625267 

-.0000014403 

6 

21.2116366294 

21.21  16359160 

—.0000007 1 32 

7 

24.3524715304 

24.3524711374 

-.0000003927 

8 

27.4934791317 

27. 4934788980 

-.0000002336 

9 

30.6346064676 

30.6346063204 

-.0000001470 

10 

33.7758202134 

33.7758201 146 

-.0000000984 

1  1 

36.9170983527 

36. 9170982856 

-.0000000665 

12 

40.0584257636 

40.0584257152 

-.0000000478 

13 

43.1997917127 

43. 1997916764 

-.0000000355 

14 

46.3411883712 

46. 3411883442 

-.0000000261 

15 

49.4826098960 

49.4826098764 

-.0000000190 

16 

52.6240518400 

52.6240518242 

-.0000000158 

17 

55.7655107528 

55.7655107416 

-.0000000110 

18 

58.9069839241 

58.9^69839139 

-.0000000099 

19 

62.0484691886 

62. 0484691802 

-.0000000082 

20 

65.1899647992 

65. 1899647918 

-.0000000070 

I 
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TABLE  VII 

CUMULATIVE  ERRORS  FOR  'x  ***  Jq(X) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  OROER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 

AFTER  THE  FIRST  ROOT  IS  COMPUTED,  THE  NEXT  ROOT  IS  COMPUTED 


FROM 

FROM 

THE  VALUE  OF  X 

THE  FIRST  ROOT. 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO* 

ROCT 

ROOT 

1 

5.5200781103 

5. 5336719233 

.0135938131 

2 

8.6537279126 

8.6700534131 

*0163255008 

3 

11.7915344390 

1 1. 8088337616 

.01 72993229 

4 

14.9309177082 

14.9486705526 

.0177528444 

5 

18.0710639674 

18. 0390637 ‘50 

*01 79997976 

6 

21.2116366294 

21.2297854787 

.0181483496 

7 

24.3524715304 

24.  3707171553 

.0182456251 

8 

27.4934791317 

27.5117911105 

.0183119791 

9 

30.6346064676 

30. 6529659065 

.0183594390 

10 

33.7753202134 

33.7942147609 

.0183945483 

11 

36.9170983527 

36. 9355196003 

.0184212479 

12 

40.0584257636 

40.0768677853 

.0184420225 

13 

43.1997917127 

43.2182502151 

.0184585024 

14 

46*3411063712 

46.3596601672 

.0184717962 

15 

49.4826098960 

49. 5010925708 

.0184826752 

16 

52.6240518400 

52.6425435262 

.0184916886 

17 

55.7655107528 

55. 7840099940 

.0184992421 

18 

58.9069839241 

58.9254095560 

.0185056320 

19 

62.0484691886 

62.0669802753 

.0185110874 

20 

65.1899647992 

65.2084805798 

.0185157806 
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TABLE  VIII 

CUMULATIVE  ERRORS  FOR  X  ***  JQ(X) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


AFTER 

THE  FIRST  ROOT 

IS  COMPUTEO,  THE  NEXT 

ROOT  IS  COMPUTED 

FROM 

THE  VALUE  OF  X 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM  THE  FIRST  ROOT. 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

5.5200781103 

5.5217120022 

.0016338920 

2 

8.6537279126 

8. 6555185516 

.0017906391 

3 

1 1 .7915344390 

1  1.7933614515 

.0018270126 

4 

14.9309177082 

14.9327572957 

.0018395875 

5 

18.0710639674 

18. 0729090152 

.0018450481 

6 

21.2116366294 

21.2134844260 

.0018477967 

7 

24.3524715304 

24. 3543208619 

.0018493315 

8 

27.4934791317 

27.4953293884 

.0018502570 

9 

30.6346064676 

30.6364573166 

.0018508493 

10 

33.7758202134 

33.7776714582 

.0018512453 

11 

36.9170983527 

36. 9189498741 

.0018515220 

12 

40.0584257636 

40. 0602774834 

.0018517202 

13 

43.1997917127 

43.2016435778 

.0018518660 

14 

46.341 1883712 

46.3430403462 

.0018519759 

15 

49.4826098960 

49.4844619567 

.0018520612 

16 

52.6240518400 

52. 6259039668 

.0018521272 

17 

55.7655107528 

55.7673629336 

.0018521811 

18 

58.9069839241 

58. 9088361468 

.0018522234 

19 

62.0484691886 

62.0503214458 

.0018522573 

20 

65.1899647992 

65.  1918170843 

.0018522854 

TABLE  IX 

CUMULATIVE  ERRORS  FOR  X  *«*  JqIX) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 

AFTER  THE  FIRST  ROOT  IS  COMPUTED,  THE  NEXT  ROOT  IS  COMPUTEO 
FROM  THE  VALUE  CF  X  LAST  OBTAINED.  ERROR  IS  CUMULATIVE 
FROM  THE  FIRST  ROOT. 


ROOT 

NO. 

TABULATED 

ROCT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5. 5172325156 

-.0028455946 

2 

8.6537279126 

8.6506294168 

-.0030984957 

3 

11.7915344390 

1  1.7883782729 

-.0031561660 

4 

14.9309177082 

14.9277417252 

-.0031759829 

5 

18.0710639674 

18.0678793993 

-.0031845679 

6 

21.2116366294 

21.2084477427 

-.0031888865 

7 

24.3524715304 

24.3492802321 

-.0031912980 

8 

27.4934791317 

27.4902863791 

-.0031927526 

9 

30.6346064676 

30. 6314127846 

-.0031936828 

10 

33.7758202134 

33.7726259045 

-.0031943080 

11 

36.9170983527 

36.9139036089 

-.0031947432 

12 

40.0584257636 

40.0552307060 

-.0031950568 

^  13 

43.1997917127 

43.  1965964232 

-.0031952894 

14 

46.3411883712 

46.3379929066 

-.0031954644 

15 

49.4826098960 

49. 4794142973 

-.0031955985 

16  . 

52.6240518400 

52.6208561352 

-.0031957048 

17 

55.7655107528 

55.7623149641 

-.0031957879 

18 

58.9069839241 

58.9037880665 

-.0031958568 

19 

62.0484691886 

62.0452732751 

-.0031959134 

20 

65.1899647992 

65.  1867688391 

-.0031959599 

I 
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TABLE  X 

CUMULATIVE  ERRORS  FOR  XI  *#*  Jq(X) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ZERO  FOR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


AFTER  THE  FIRST  ROOT 

IS  COMPUTED,  THE  NEXT 

ROOT  IS  COMPUTED 

FROM 

THE  VALUE  OF  XI 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM 

THE  FIRST  ROOT. 

ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5254224626 

.0053443524 

2 

8.6537279126 

8.6603153765 

.0065874640 

3 

11.7915344390 

1  1. 7985812158 

.0070467769 

4 

14.9309177082 

14.9381818194 

.0072641112 

5 

18.0710639674 

1 8. 078447483S 

.0073835164 

6  • 

21.2116366294 

21.  2190926266 

.0074559972 

7 

24.3524715304 

24.3599747745 

.0075032441 

8 

27.4934791317 

27.5010148650 

.0075357334 

9 

30.6346064676 

30. 6421654909 

.0075590234 

10 

33.7758202134 

33.7833964955 

.0075762823 

1  1 

36.9170983527 

36.9246877786 

.0075894268 

12 

40.0584257636 

40.  0660254285’ 

.0075996651 

13 

43.1997917127 

43.2073995080 

.0076077954 

14 

46.3411883712 

46. 3488027304 

.0076143594 

15 

49.4826098960 

49.4902296308 

.0076197354 

16 

52.6240518400 

52.6316760313 

.0076241914 

17 

55.7655107528 

55. 7731386814 

.0076279288 

18 

58.9069839241 

58.9146150155 

.0076310916 

19 

62.0484691886 

62.0561029809 

.0076337926 

20 

65.1899647992 

65.  19760091 60 

.0076361171 
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TABLE  XI 

CUMULATIVE  ERRORS  FOR  71  ***  Jq ( X ) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  OROER  ZERO  FpR  ROOTS 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


AFTER  THE  FIRST  ROOT 

IS  COMPUTED,  THE  NEXT 

ROOT  IS  COMPUTED 

FROM 

THE  VALUE  OF  X2 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM 

THE  FIRST  ROOT. 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

5.5200781103 

5.5194687883 

-.0006093219 

2 

8.6537279126 

8.6530707316 

-.0006571810 

3 

11.7915344390 

11. 7908666495 

-.0006677894 

4 

14.9309177082 

14.9302463082 

-.0006714000 

5 

18.0710639674 

18.0703910086 

-.0006729586 

6 

21.2116366294 

21. 2109628869 

-.0006737421 

7 

24.3524715304 

24.3517973502 

-.0006741801 

8 

27.4934791317 

27.4928046870 

-.0006744443 

9 

30.6346064676 

30.6339318538 

-.0006746134 

10 

33.7758202134 

33.7751454851 

-.0006747281 

1  1 

36.9170983527 

36.9164235452 

-.0006748069 

12 

40.0584257636 

40. 0577508984 

-.0006748646 

13 

43.1997917127 

43. 1991168046 

-.0006749080 

14 

46.3411883712 

46. 3405134305 

-.0006749401 

15 

49.4826098960 

49.481934931 1 

-.0006749644 

16 

'  52.6240518400 

52.6233768547 

-.0006749846 

17 

55.7655107528 

55.7648357525 

-.0006749997 

10 

58.9069839241 

58.9063089108 

-.0006750130 

19 

62.0484691686 

62.0477941642 

-.0006750239 

20 

65.1899647992 

65.1892897654 

-.0006750335 
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TABLE  XII 

VALUES  OF  X*  ***  Jj (X) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE _ 

p(xim  =  i+n-4N2)/4(x(i  ))2  xu  +  i)  =  a+H/  i/p(x(m 


ERROR 

IS  THE  DIFFERENCE 

8ETWEEN  TABULATED 

AND  COMPUTED  ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

6.9976376109 

-.0179490588 

2 

10.1734681348 

10. 1686351451 

-.0048329895 

3 

13.3236919360 

13.3217202593 

-.0019716767 

4 

16.4706300506 

16. 4696368496 

-.0009932009 

5 

19.6158585101 

19.6152890925 

-.0005694173 

6 

22.7600853802 

22. 7597279320 

-.0003574477 

7 

25.9036720875 

25.9034342724 

-.0002378150 

8 

29.0468205335 

29.0466620042 

-.0001665289 

9 

32.1896799095 

32.  1895587798 

-.0001211293 

10 

35.3323075492 

35. 3322166996 

r. 0000908489 

11 

38.4747662339 

38.4746963540 

-.0000698795 

12 

41.6170942122 

41. 6170393117 

-.0000549003 

13 

44.7593189972 

44.7592750816 

-.0000439151 

14 

47.9014608869 

47.9014252098 

-.0000356767 

15 

51.0435351832 

51.0435058065 

-.0000293767 

16 

54.1855536401 

54. 1855291631 

-.0000244770 

17 

57.3275254359 

57.3275048267 

-.0000206090 

18 

60.4694578443 

60. 4694403261 

-.0000175173 

19 

63.6113566970 

63.61 13416851 

-.0000150119 

20 

66.7532267310 

66.7532137688 

-.0000129620 

44 

I 


TABLE  XIII 

VALUES  OF  9^  ***  JjU) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


CHI  IS  THE  VALUE  X  MUST  ASSUME  TO  CONVERGE  TO  B 
{  {  8- A  ) /2)\y  (  1-4N2)/{  {  B-A  )2-7l2)  ) 


ERROR 

IS  THE  DIFFERENCE 

BETWEEN  X  AND  CHI 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROCT 

ROOT 

1 

7.0155866698 

5.3312873839 

1.6663502270 

2 

10.1734681348 

8. 5374782400 

1.6311569051 

3 

13.3236919360 

1  1.7071  112888 

1.6146089705 

4 

16.4706300506 

14.  8645464401 

1.6050904095 

5 

19.6158585101 

18.0163550223 

1.5989340702 

6 

22.7600853802 

21.  1610798943 

1.5986480378 

7 

25.9036720875 

24. 3119656080 

1.5914686644 

8 

29.0463285335 

27. 4576333771 

1.5890286271 

9 

32.1896799095 

30.6024408266 

1.5871179532 

10 

35.3323075492 

33. 7466444736 

1.5855722260 

11 

38.4747662339 

36. 8904200587 

1.5842762953 

12 

41.6170942122 

40. 0338375270 

1.5832017846 

13 

44.7593189972 

43.  1769980257 

1.5822770558 

14 

47.9014608869 

46. 3199358247 

1.5814893851 

15 

51.0435351832 

49.4627229534 

1.5807828531 

16 

54.1 B5553640 1 

52.  6053818800 

1.5801472832 

17 

57.3275254359 

55.7479370888 

1.5795677379 

18 

60.4694578443 

58.  8902029432 

1.5792373829 

19 

63.6113566970 

62.0326487795 

1.5786929056 

20 

66.7532267310 

65.  1751772780 

1.5780364908 
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TABLE  XIV 
VALUES  OF  X  *** 


J  ^  (  X  ) 


ZEROS  OF  THE  BESSEL  FUNCTION 
ONE  THROUGH  TWENTY,  BASED  ON 
P(X(I)J  =  1+!1-4N2)/4I(A+X(I) 
ERROR  IS  THE  DIFFERENCE  BETWE 


OF  ORDER  ONE  FOR  ROOTS  NUMBER 
THE  RECURSION  FORMULAE 
>/2)2  7(1+1)  =  A+77/  /PCX!  I  )  J 

EN  TABULATED  AND  COMPUTED  ROOTS 


ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7.0141418219 

-.0014448478 

2 

10.1734681348 

10. 1732513588 

-.0002167758 

3 

13.3236919360 

13. 3236309085 

-.0000610273 

4 

16.4706300506 

1  6.  4706066190 

-.0000234314 

5 

19.6158585101 

19.6158476649 

-.0000108451 

6 

22.7600853802 

22. 7600786900 

-.0000066901 

7 

25.903672087  5 

25. 9036688204 

-.0000032669 

8 

29.0468285335 

29.0468265270 

-.0000020061 

9 

32.1896799095 

32.  1896786084 

-.0000013006 

10 

35.3323075492 

35.3323066691 

-.0000008800 

1  1 

38.4747662339 

38.4747656174 

-.0000006158 

12 

41.6170942122 

41. 6170937680 

-.0000004442 

13 

44.7593189972 

44. 7593186684 

-.0000003282 

14 

47.9014608869 

47. 9014606383 

-.0000002479 

15 

51.0435351832 

51. 0435349923 

-.0000001903 

16 

54.1855536401 

54. 1855534911 

-.0000001483 

17 

57.3275254359 

57.  3275253186 

-.0000001170 

18 

60.4694578443 

60. 4694577483 

-.0000000957 

19 

63.6113566970 

63.6113566197 

-.0000000766 

20 

66.7532267310 

66.7532266676 

-.0000000615 
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TABLE  XV 

VALUES  OF  X  ***  J  ^  (  X ) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 

ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE _ 

=  1  +  (1  -  4N2)/  4AXII)  $(I+1)  *  A  +7T/\/p  dim  ) 
ERROR  IS  THE  DIFFERENCE  BETWEEN  TABULATED  AND  COMPUTED  ROOTS 


ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7. 0180468041 

.0024601344 

2 

10.1734681348 

10. 1738162260 

.0003480912 

3 

13.3236919360 

13.3237883153 

.0000963795 

U 

16.4706300506 

16.4706667955 

.0000367450 

5 

1  9. 6 1  5P  53  51  01 

19.6158754556 

.0000169458 

6 

22.7600853802 

22.7600932522 

.0000078722 

7 

25.9036720075 

25.9036771739 

.0000050868 

8 

29.0468285335 

29.0468316572 

.0000031241 

9 

32.1 S96799095 

32. 1896819314 

.0000020222 

10 

35.3323075492 

35.3323089145 

.0000013660 

11 

38.4747662339 

38.4747671895 

.0000009564 

12 

41.6170942122 

41. 6170949005 

.0000006892 

13 

44.7593189972 

44.7593195057 

.0000005093 

14 

47.9014608869 

47.9014612706 

.0000003842 

15 

51.0433351832 

51.0435354784 

.0000002958 

16 

54.1055536401 

54. 1855538711 

.0000002316 

17 

57.3275254359 

57.  3275256194 

.0000001642 

18 

60.4694578443 

60.4694579896 

.0000001462 

19 

63.6113566970 

63.6113568163 

.0000001198 

20 

66.7532267310 

•66.  7532.268297 

.0000000997 

TABLE  XVI 

VALUES  OF  Xl  ***  J^(X) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FC)R  ROOTS  NUMBER 

ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE _ 

P(F»(I))  ”  (  P(  X(  I  )  )  +  P  (X(  I  )  )  )/ 2  xT(l)  *  A+7C/V  PCxTt  I  )  ) 


ERROR 

IS  THE  DIFFERENCE 

BETWEEN  TABULATED 

AND  COMPUTED  ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7. 0077901660 

-.0077965038 

2 

10.1  73468  1 3  4  ti 

10. 17 1221693C 

-.0022464416 

3 

13.3236919360 

13. 3227535994 

-.0009383365 

4 

1 6.470630C506 

16.4701516442 

-.0004784064 

5 

19.6158585101 

19. 6155822147 

-.0002762953 

.6 

22.7600853802 

22.7599105686 

-.00017481 14 

7 

25.9036720875 

25.9035557127 

-.0001163747 

8 

29.0468285335 

29.0467468258 

-.0000817076 

9 

32.1896799095 

32. 1896203524 

-.0000595564 

10 

35.3323075492 

35.3322628057 

-.0000447430 

1 1 

38.4747662339 

38.4747317713 

-.0000344625 

12 

41.6170942122 

41.6170671051 

-.0000271062 

13 

44.7593189972 

44. 7592972936 

-.0000217033 

14 

47.9014608869 

47.9014432402 

-.0000176465 

15 

51.0435351832 

51.0435206424 

-.0000145406 

16 

54.1855536401 

54. 1855415171 

-.0000121228 

17 

57.3275254359 

57. 3275152231 

-.0000102125 

18 

60.4694578443 

60.4694491578 

-.0000086856 

19 

63.6113566970 

63. 6113492502 

-.0000074461 

20 

66.7532267310 

66.7532202993 

-.0000064312 

TABLE  XVII 

VALUES  OF  X2  «**  JjIX) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 

ONE  THROUGH  TWENTY,  BASEO  ON  THE  RECURSION  FORMULAE _ 

pixFim  =  (p(xm)+p($m))/2  x2(i>  =  a+F/ vp i Yz m  ) 


ERROR 

IS  THE  DIFFERENCE 

BETWEEN  TABULATED 

AND  COMPUTED  ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7. 0160913781 

.0005047083 

2 

10.1734681348 

10. 1735337211 

.0000655864 

3 

13.3236919360 

13. 3237096060 

.0000176702 

4 

16.4706300506 

16. 4706367063 

.0000066559 

5 

19.6158585101 

19.6158615602 

.0000030501 

6 

22.7600853802 

22.  76008597  11 

.0000005910 

7 

25.9036720875 

25. 9036729969 

,0000009099 

8 

29.0468285335 

29. 0468290923 

.0000005590 

9 

32.1896799095 

32.  1896802699 

.0000003608 

10 

35.3323075492 

35. 3323077913 

.0000002430 

1  1 

38.4747662339 

38.4747664034 

.0000001703 

12 

41.6170942122 

41. 6170943342 

.0000001225 

13 

44.7593189972 

44.7593190875 

.0000000905 

.  14 

47.9014608869 

47. 9014609549 

.0000000681 

15 

51.0435351832 

51.0435352353 

.0000000527 

16 

54.1855536401 

54.  185553681  1 

.0000000416 

17 

57.3275254359 

57.3275254695 

.0000000336 

18 

60.4694578443 

60. 4694578694 

.0000000253 

19 

63.6113566970 

63.6113567185 

.0000000216 

20 

66.7532267310 

66. 7532267496 

.0000000191 

TABLE  XVIII 

CUMULATIVE  ERRORS  FOR  'J?  ***  J^X) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 

AFTER  THE  FIRST  ROOT  IS  COMPUTED,  THE  NEXT  ROOT  IS  COMPUTED 


FROM 

THE  VALUE  OF  ^ 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM 

THE  FIRST  RCOT . 

ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

I 

7.0155866698 

6. 9976376109 

“.01 79490588 

2 

10.1734681348 

10. 1507267661 

-.0227413687 

3 

13.3236919360 

13.2990017352 

-.0246902007 

4 

16.4706300506 

16. 4449597476 

-.0256703027 

5 

19.6158585101 

19.5896266408 

-.0262316693 

6 

22.7600853802 

22. 7335015247 

- . 0265838551 

7 

25.9036720875 

25. 8768550311 

-.0268170560 

8 

29.0468285335 

29.  0198475332 

-.0269810002 

9 

32.1896799095 

32. 1625796901 

-.0271002193 

10 

35.3323075492 

35. 3051179312 

-.0271896177 

11 

38.4747662339 

38. 4475078629 

-.0272583706 

12 

41.6170942122 

41. 5897818329 

-.0273123789 

.  13 

44.7593189972 

44. 7319634212 

-.0273555757 

14 

'47.9014608869 

47. 8740702206 

-.0273906654 

15 

51.0435351832 

51.  0161  156254 

-.0274195571 

16 

54.1855536401 

54.  1581  1001  14 

-.0274436282 

17 

57.32  /'S 254359 

57.  3000615416 

-.0274638942 

18 

60.4694578443 

60.4419767251 

-.0274811186 

19 

63.61 13566970 

63.  5838608174 

-.0274958793 

20 

66.7532267310 

66.7257181071 

-.0275086237 
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TABLE  XIX 

CUMULATIVE  ERRORS  FOR  X  **»  J,{X) 

^  \ 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 

AFTER  THE  FIRST  ROOT  IS  COMPUTED,  THE  NEXT  ROOT  IS  COMPUTED 


FROM 

THE  VALUE  OF  X 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM 

THE  FIRST  R COT • 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7.0141418219 

-.0014448478 

2 

10.1734681348 

10. 1718119476 

-.0016561872 

3 

13.3236919360 

13.  3219771462 

-.00 17147896 

4 

1 6. 470630C506 

16.4688930572 

-.0017369931 

5 

19.6158585101 

19.6141113704 

-.0017471395 

6 

22.7600853802 

22.7583319838 

-.0017533960 

7 

25.9036720875 

25. 9019167111 

-.0017553763 

8 

29.0468285335 

29. 0450713504 

-.0017571828 

9 

32.1896799095 

32. 1879215697 

-.0017583393 

10 

35.3323075492 

35.3305484364 

-.0017591120 

1  1 

38.4747662339 

38. 4730065875 

-.0017596461 

12 

41.6170942122 

41. 6153341858 • 

-.00  17600260 

13 

44.7593189972 

44.7575586932 

-.0017603031 

14 

47.9014608869 

47. 8997003762 

-.0017605102 

15 

51.0435351832 

51. 04177451 60 

-.0017606668 

16 

54. 1855536401 

54. 1837928528 

-.00 17607870 

17 

57.32/5254359 

57. 3257645555 

-.0017608804 

18 

60.4694578443 

60.4676968884 

-.0017609558 

19 

63.611 3566970 

63. 6095956815 

-.0017610151 

20 

66.7532267310 

66.7514656689 

-.0017610620 

TABLE  XX 

CUMULATIVE  ERRORS  FOR  £  ***  JjIX) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


AFTER 

THE  FIRST  ROOT 

IS  COMPUTEO,  THE  NEXT 

ROOT  IS  COMPUTED 

FROM 

THE  VALUE  OF  $ 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM  THE  FIRST  ROOT. 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7.0180468041 

.0024601344 

2 

10. 1734681348 

10.  1762664455 

.0027983109 

3 

13.3236919360 

13. 3265823792 

.0028904433 

4 

16.4706300506 

16. 4735551225 

.0029250719 

5 

19.6158585101 

19.6187993325 

.0029408227 

6 

22.7600853802 

22. 763033336! 

.0029479564 

7 

25.9036720875 

25. 9066256429 

.0029535557 

8 

29.0468285335 

29. 0497848745 

.0029563413 

9 

32. 1 896799095 

32.  1926380284 

.0029581190 

TO 

35.3323075492 

35. 3352668518 

.0029593028 

1  1 

38.4747662339 

38.4777263533 

.0029601197 

12 

41.6170942122 

41. 6200549109 

.0029606996 

13 

44.7593189972 

44.7622801 177 

.0029611211 

14 

47.9014608869 

47. 9044223214 

.0029614345 

15 

51.043535 1832 

51. 0464968551 

.0029616724 

16 

54.1855536401 

54. 1885154955 

.0029618555 

17 

57.3275254359 

57.3304874348 

.0029619992 

16 

60.4694578443 

60. 4724199548 

.0029621108 

19 

63.6113566970 

63. 6143188979 

.002962201 1 

20 

66.7532267310 

66. 7561890054 

.0029622754 

) 
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TA8LE  XXI 

CUMULATIVE  ERRORS  FOR  XI  ***  J^JX) 

ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FQR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 

AFTER  THE  FIRST  ROOT  IS  COMPUTED,  THE  NEXT  ROOT  IS  COMPUTEO 


FROM 

FROM 

THE  VALUE  OF  XI 

THE  FIRST  ROOT. 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7.0077901660 

-.0077965038 

2 

10.1734681348 

10. 1634497307 

“.01 00184040 

3 

13.3236919360 

13.3127478317 

-.0109441041 

4 

16.4706300506 

16.4592144503 

-.01 14155999 

5 

1 9.6158585101 

19.6041707369 

-.01 16877730 

6 

22.7600853802 

22.7482254347 

-.01 18599452 

T 

25.9036720875 

25.8916985528 

-.01 19735343 

8 

29.0468285335 

29.0347745535 

-.0120539797 

9 

32.1896799095 

32.  1775672976 

-.0121126117 

10 

35.3323075492 

35.  3201508904 

-.0121566580 

1 1 

38.4747662339 

38.4625756508 

-.0121905831 

12 

41.6170942122 

41. 6048769467 

-.0122172653 

13 

44.7593189972 

44.7470803680 

-.0122386289 

14 

47.9014608869 

47.  8892048877 

-.0122559989 

15 

51.0435351832 

51.0312648714 

-.0122703116 

16 

54.1855536401 

54. 1732713953 

-.0122822442 

1? 

57.3275254359 

57.3152331384 

-.0122922968 

18 

60.4694578443 

60.4571569972 

-.0123008466 

19 

63.6113566970 

63.5990485204 

-.0123081761 

20 

66.7532267310 

66.7409122232 

-.0123145066 

/ 
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TABLE  XXII 

CUMULATIVE  ERRORS  FOR  X2  ***  JjIX) 


ZEROS  OF  THE  BESSEL  FUNCTION  OF  ORDER  ONE  FOR  ROOTS  NUMBER 
ONE  THROUGH  TWENTY,  BASED  ON  THE  RECURSION  FORMULAE 


AFTER  THE  FIRST  ROOT 

IS  COMPUTEO,  THE  NEXT 

ROOT  IS  COMPUTED 

FROM 

THE  VALUE  OF  X2 

LAST  OBTAINED.  ERROR 

IS  CUMULATIVE 

FROM 

THE  FIRST  ROOT. 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7. 0160913781 

.0005047083 

2 

10.1734681348 

10. 1740364628 

.0005683281 

3 

13.3236919360 

13. 3242770866 

.0005851508 

4 

1 6. 470630  0506 

16.4712214330 

.0005913828 

5 

19.6158585101 

1 9. 6 1 6452702c 

.0005941927 

6 

22. 7600853802 

22. 76G6000 1 52 

.0005946353 

7 

25.9036720875 

25.9042685339 

.0005964467 

8 

29. 0468285335 

29. 0474254708 

.0005969373 

9 

32.1896799095 

32. 1902771583 

.0005972489 

10 

35.3323075492 

35.  3329050038 

.0005974551 

1  1 

38.4747662339 

38.4753638301 

.0005975968 

12 

41.6170942122 

41. 6176919080 

.0005976967 

13 

44.7593189972 

44.7599167656 

.0005977688 

14 

47. 9  u. *.606  86  9 

47. 9020587085 

.0005978224 

15 

51. 043535 1 832 

51.0441330457 

.0005978626 

16 

54.1855536401 

54. 1861515343 

.0005978945 

17 

57.0275254359 

57.3281253553 

.0005979196 

18 

60.4694578443 

60.4700557813 

.0005979377 

19 

63.6113566970 

63.61  19546490 

.0005979527 

20 

66.7532267310 

66.7538246959 

.0005979659 

i 
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